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Conferences 

Progress  on  subject  contract.  Watkins-Johnson 

Company  personnel  and  Mr.  Robert  M.  True  of 
ERADCOM,  Fort  Monmouth,  N.J.  Held  at 
Watkins-Johnson  Company  on  20-22  March  1979. 


Program  for  the  Next  Internal 

The  Program  Plan  shown  in  Figure  2-2  and  described  in  Section  6.0  represents 
our  best  estimate  of  work  to  be  carried  out  during  the  next  reporting  period. 
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1. 1 Objective 

The  program  objective  is  to  reduce  the  output  power  level  of  spurious 
noise  signals,  intermodulation  (IM)  products  and  harmonic  distortion 
generated  by  deflection  modulated  electron  beam  semiconductor  (EBS) 
amplifiers. 

1.2  Technical  Approach 

Figure  1-1  illustrates  the  configuration  of  a deflected  beam  EBS  amplifier. 
This  type  of  amplifier  has  been  developed  by  Watkins-Johnson  Company 
over  a period  of  several  years.  Measurements  of  existing  EBS  ampli- 
fiers will  be  made  to  determine  the  typical  values  of  IM  products , spuri- 
ous noise  and  harmonic  distortion.  The  existing  electron  beam  profile 
will  be  characterized  using  a slit  beam  analyzer.  The  EBS  performance 
will  be  correlated  to  the  measured  beam  profile  using  a i.'.athematical 
analysis  implemented  by  a computer  program.  A second  computer  pro- 
gram will  perform  an  analysis  of  the  expected  beam  profile  generated  by 
the  existing  gun  geometry  and  these  results  will  be  correlated  with  the 
measured  beam  profile.  Modifications  will  then  be  performed  on  the 
electron  gun  to  improve  the  linearity,  IM  products  and  harmonic  distor- 
tion; and  the  re-designed  gun  will  be  fabricated  and  tested  on  the  beam 
analyzer.  Two  devices  using  the  re-designed  gun  will  be  fabricated, 
tested  and  delivered  as  part  of  this  contract. 
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Figure  1-1.  Deflected  Beam  EBS  Configuration 
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2.0  PROGRAM  SUMMARY 

2. 1 Summary  of  Work 

During  the  four  month  period  covered  by  this  report , the  following  was 
accomplished : 

1.  The  existing  EBS  amplifiers  having  2x6  diode  arrays  were 
completely  characterized  in  terms  of  linearity  and  intermodulation 
products  as  a function  of  device  efficiency. 

2.  The  computer  aided  analysis  of  the  present  gun  performance  was 
completed.  Several  conclusions  were  drawn  from  the  results,  and 
areas  of  potential  improvements  were  identified. 

3.  Work  continued  on  methods  of  analyzing  actual  shape  of  electron 
beams.  Calibration  of  an  existing  slit  diode  array  was  not  satis- 
factory, so  a fabrication  of  a new  slit  diode  array  was  begun. 
Phosphor  screens  were  used  as  a secondary  method  of  beamshape 
observation . 

4.  Several  guns  were  checked  using  the  phosphor  screen  method. 
Poor  beam  shape  was  observed,  and  the  guns  were  checked  for 
dimensional  conformation  to  the  present  design.  Some  problems 
were  identified  in  this  regard,  and  redesigned  tooling  and  parts 
were  ordered  to  correct  these  problems. 

2.2  Program  Schedule 

Figure  2-1  shows  the  Program  Schedule,  updated  as  of  April  1979.  The 
primary  cause  of  slippage  has  been  lack  of  a suitable  method  of  measur- 
ing beam  shape.  This  problem  has  been  solved  by  the  use  of  phosphor 
screens  and  the  fabrication  of  a new  slit  diode  array. 

All  other  tasks  are  proceeding  on  schedule  with  one  exception:  because 
of  results  obtained  by  beam  shape  analysis  performed  as  part  of  this 
program,  emphasis  has  been  shifted  away  from  the  development  of  the 
metal-ceramic  version  of  the  sheet  beam  gun  and  increased  emphasis  has 
been  put  on  fabrication  techniques  using  the  present  glass  envelope  con- 
figuration. For  this  reason,  the  two  tasks  having  to  do  with  develop- 
ment of  the  metal  ceramic  gun  have  been  deleted  from  the  program. 
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Figure  2-1.  Program  Milestone  Schedule 
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CHARACTERIZATION  OF  EXISTING  DEVICES 


Extensive  linearity  and  intermodulation  measurements  were  performed  on 
a deflected  beam  amplifier,  WJ- 3662-1  S/N  9.  The  device  was  set  up  for 
saturated  output  power  of  46.0  dBm  at  maximum  efficiency.  The  balanced 
two  tone  method  of  measurement  was  used,  and  in  most  cases,  the  funda- 
mental tones  and  their  3rd,  5th,  and  7th  intermodulation  products  were 
measured  over  a 20-33  dB  dynamic  range.  The  parameter  varied  during 
these  tests  was  the  voltage  on  anode  3,  the  next  to  last  anode  in  the  gun 
stack.  (Anode  4 is  held  at  ground  potential.)  Figures  3-1  through  3-5 
summarize  these  results,  and  Figure  3-6  expresses  the  tradeoff  between 
linearity  and  efficiency.  The  criterion  for  linearity  is  the  maximum  3rd 
order  intermodulation  product  referenced  to  fundamental  saturation  level; 
note  that  the  case  for  "best"  linearity  by  this  criterion  (A3  = 3400)  is 
not  coincident  with  the  case  where  the  fundamental  deviates  the  least 
from  constant  gain  (A3  = 3500),  nor  with  the  case  where  single  tone 
saturated  efficiency  is  maximum  (A3  = 3300).  The  reason  for  this  is 
attributed  to  the  sidelobes  present  in  the  beam  density  profile,  which  has 
been  observed  on  phosphor  screens.  At  maximum  efficiency,  the  amplifier 
is  running  in  class  C,  i.e.  , with  a tightly  focused  beam  smaller  than  the 
diode  spacing.  At  other  points,  class  AB  predominates,  with  current 
spillover  severely  reducing  efficiency.  Since  the  beam  profile  as  pre- 
dicted by  computer  for  class  C is  not  perfectly  rectangular,  the  inter- 
modulation products  are  more  severe  than  for  class  AB  , which  in  spite 
of  spillover  has  a better  beamshape  within  the  width  of  the  diode  spacing. 
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Figure  3-2.  IMD  Performance  (A„  = 3200  V) 


Figure  3-3.  IMD  Performance  for  the  Case  Where  Maximum  Target  Efficiency 
Is  Attained  at  Single  Tone  Saturation  (A„  = 3300  V). 


Figure  3-5.  IMD  Performance  (A„  = 3500  V). 
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4.0  MEASUREMENT  OF  ELECTRON  BEAM  SHAPE 

4. 1 Slit  Diode  Array 

EBS  device  linearity  is  largely  dependent  upon  the  profile  of  the  current 
density  of  the  sheet  electron  beam.  An  accurate  means  is  required  to 
measure  this  profile.  A special  diode  target  called  a "slit  beam  analyzer" 
shown  schematically  in  Figure  4-1  was  prepared  for  this  purpose.  The 
target  contains  10  diodes,  each  having  a metalized  top  surface  with  a 
narrow  slit  etched  through  the  metalization  to  the  active  diode  layer.  If 
an  electron  beam  is  swept  across  this  array  as  illustrated , each  diode 
will  respond  only  to  the  portion  of  the  beam  striking  the  narrow  slit  area. 
The  diode  current  from  each  of  the  diodes  in  the  array  can  be  monitored 
and  a profile  of  the  beam  density  can  be  thus  obtained. 

During  the  first  four  months  of  the  contract , an  existing  slit  diode  array 
was  assembled  to  a sheet  beam  electron  gun  and  was  tested. 

Because  of  the  inconsistencies  of  response  from  diode  to  diode,  it  was 
suspected  that  the  array  itself  was  not  responding  uniformly  to  excita- 
tion by  the  electron  beam.  It  became  apparent  that  a calibration  of  the 
slit  beam  analyzer  array  was  required  before  meaningful  data  could  be 
obtained. 

A gun  having  a well-defined  circular  beam  of  small  radius  was  needed  to 
calibrate  the  slit  diode  array.  Such  a gun  is  typically  used  in  a cathode 
ray  tube  (CRT)  for  information  display.  The  Stewart  Division  of  Watkins- 
Johnson  Company  produces  these  "pencil  beam"  guns  for  CRT  use,  and 
arrangements  were  made  to  purchase  such  a gun  from  them.  A suitable 
production  type  gun  was  identified  and  the  problems  associated  with 
interfacing  the  gun  to  the  slit  beam  analyzer  target  were  solved.  In 
order  to  provide  sufficient  deflection  of  the  electron  beam  to  cover  the 
entire  target,  using  available  sawtooth  generators,  the  deflection  coil  to 
be  used  must  be  designed  to  match  the  physical  dimensions  of  the  gun 
and  its  coil  impedance  should  be  roughly  the  same  as  the  sawtooth  gene- 
rator output  impedance.  A suitable  coil  was  designed  and  procured  from 
an  outside  vendor. 

When  the  pencil  beam  gun  was  operated  with  the  slit  diode  array,  the 
diode  current  obtained  was  insufficient  to  provide  a high  enough  signal- 
to- noise  ratio.  In  addition,  problems  were  experienced  in  obtaining  a 
reasonable  amount  of  cathode  current  from  the  gun. 

The  gun  was  removed  from  the  slit  diode  array,  and  sent  back  to  the 
Stewart  Division  for  installation  of  a new  cathode.  In  conjunction  with 
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no  beam  deflection. 


•Direction  of  Beam  Deflection 


Figure  4-1.  Slit  Beam  Analyzer  Array 
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this,  fabrication  of  a new  slit  diode  array  was  begun,  since  several  of 
the  diodes  in  the  slit  diode  array  were  not  functioning. 

The  new  slit  diode  array  will  have  wider  (0.002  inches)  slits  to  provide 
greater  response  amplitude.  Since  the  optimum  beam  width  is  0.040 
inches,  resolution  will  still  be  far  in  excess  of  what  is  required.  The 
new  slit  diode  array  will  be  first  tested  with  the  pencil  beam  gun  to  cali- 
brate it,  and  can  then  be  used  to  test  beam  shape  on  EBS  sheet  beam 
guns. 

. 2 Phosphor  Screen  Beam  Measurements 

Because  of  the  problems  experienced  with  the  slit  diode  array,  an  alter- 
nate method  of  testing  the  beam  shape  was  explored.  A phosphor  screen 
similar  to  that  used  on  a CRT  is  mounted  to  the  EBS  device  in  place  of 
the  EBS  diode  array.  The  beam  cross-section  is  observed  when  the  beam 
strikes  and  illuminates  the  phosphor.  The  phosphor  is  very  sensitive  to 
electron  bombardment  so  that  very  low  beam  current  is  sufficient  to 
observe  the  beam  shape.  However,  it  is  not  possible  to  run  high  beam 
current  at  normal  EBS  operating  levels,  because  the  phosphor  is  evapo- 
rated at  low  power  levels.  This  disadvantage  was  partially  alleviated  by 
pulsing  the  electron  beam  at  very  low  duty  cycle.  In  addition,  there 
seems  to  be  little  change  in  beam  shape  as  the  beam  current  is  increased. 

The  first  gun  tested  was  originally  used  on  WJ-3662-1  S/N  16,  which 
failed  at  low  power  levels  for  unknown  reasons.  The  beam  shape  shown 
on  the  screen  had  a much  greater  intensity  at  one  end  than  at  the  other, 
and  the  beam  was  not  straight  but  exhibited  a significant  amount  of 
curvature.  Maximum  deviation  of  the  beam  center  from  straight  line  was 
approximately  0.04  inch.  A photograph  of  the  beam  image  on  the  phosphor 
screen  is  shown  in  Figure  4-2. 

The  second  gun  tested  produced  a beam  having  better  overall  shape,  but 
still  with  too  much  curvature  and  with  some  increase  in  intensity  and  beam 
size  at  the  very  ends  of  the  beam,  shown  in  Figure  4-3.  The  third  gun 
was  the  worst  tested,  having  a deviation  from  straight  line  of  approximately 
0.05  inches,  shown  in  Figure  4-4. 

The  first  and  third  electron  guns  were  removed  from  the  phosphor  screens 
and  visually  examined.  In  both  guns,  misalignment  was  noted  in  the  plane 
of  the  beam  between  the  cathode  and  focus  electrode  assembly  and  the 
anode  assembly.  The  fixturing  used  to  assemble  the  guns  was  checked  and 
the  problem  was  apparent  - when  the  lengthwise  dimension  of  the  focus 
electrode  assembly  was  increased  in  early  1978,  the  fixture  had  not  been 
modified  to  take  up  the  slack.  New  fixturing  to  correct  the  problem  was 
ordered. 
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Figure  4.2  - Beam  shape  of  gun  #42  measured  on  a phosphor  screen 


Figure  4.3  - Beam  shape  of  gun  #47  measured  on  phosphor  screen. 
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Figure  4.4  - Beam  shape  of  gun  #49  measured  on  a phosphor  screen 
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The  last  gun,  which  had  the  most  pronounced  beam  curvature,  also  had 
a noticeable  warpage  in  the  first  anode. 

To  reduce  the  "dogbone"  effect  of  higher  beam  size  and  intensity  near 
the  ends  of  the  beams , modifications  were  begun  on  the  gun  electrode 
designs.  The  modifications  will  increase  the  length  of  the  slot  in  all  the 
gun  electrodes. 

To  increase  the  precision  of  the  anodes  themselves,  a new  method  of 
fabricating  the  anodes  was  developed.  Presently,  each  anode  is  machined 
to  the  desired  dimensions.  The  new  method  involves  machining  each 
anode  in  two  parts,  placing  the  parts  on  a mandrel,  and  laser  welding  the 
two  pieces  together.  Since  the  mandrel  can  be  made  with  very  good 
accuracy,  the  internal  dimensions  of  the  finished  anode  should  be  better 
using  the  new  method. 
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5.0  COMPUTER  AIDED  ANALYSIS  OF  ELECTRON  GUN 


During  the  present  reporting  period,  the  final  computer  program  was 
written  to  complete  the  correlation  between  the  rectangular  gun  geometry 
and  the  Fourier  components  of  the  output  signal.  All  of  the  programs 
currently  used  for  gun  analysis  are  listed  in  Table  I.  The  first  three 
perform  the  function  of  determining  beam  profile  through  the  acceleration 
region,  i.e.  , up  to  the  meanderline.  The  fourth  accounts  for  beamspread 
through  the  drift  region,  from  the  meanderline  to  the  target.  Finally, 
the  last  performs  a convolution  of  the  beam  density  profile  with  diode 
shape  to  determine  the  output  waveform,  then  does  a Fourier  transform  on 
this  waveform  to  determine  the  spectral  components  and  calculates  the  tar- 
get efficiency  for  the  fundamental. 

The  design  of  a new  gun  structure  has  typically  involved  the  use  of  only 
the  first  two  programs.  Generally,  the  quality  of  the  beamshape  is  quite 
obvious  after  the  second  program  is  run,  and  if  there  are  severe  prob- 
lems with  a design,  there  is  no  point  in  carrying  out  further  calculations. 
For  this  reason.  Parts  (1)  and  (2)  of  Table  I have  been  streamlined  to  run 
with  one  job  submission.  The  remainder  of  the  programs,  however,  take  a 
great  deal  of  effort  to  run,  and  are  used  only  to  check  out  a design  that 
has  been  finalized  from  (2).  A complete  listing  of  all  above  programs  is 
given  in  Appendices  1 through  IV . 
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6.0  MODIFICATION  TO  ELECTRON  GUN 


Figure  6-1  and  6-2  illustrate  the  2-dimensional  beam  profile  resulting 
from  the  present  electrode  geometry  and  ideal  Pierce  geometry,  respec- 
tively. Figure  6-2  has  served  as  a model  for  the  development  of  a new 
anode  geometry,  based  on  the  following  areas  of  improvement  that  have 
been  identified  from  Figure  6-1: 

1.  The  effective  cathode  location  is  not  necessarily  at  the  zero 
potential  line. 

2.  The  grid  is  biased  negative  relative  to  the  cathode,  in  order  to 
reduce  emission.  As  mentioned  previously,  this  crowds  the 
emission  density  into  the  center  of  the  cathode. 

3.  Potentials  farther  from  the  cathode  are  such  that  emission  at  the 
cathode  can  not  be  reduced  without  causing  the  beam  area  to  be 
smaller,  thereby  increasing  current  density. 

Problems  (2)  and  (3)  are  being  addressed  by  construction  of  a gridded 
anode  in  front  of  the  cathode,  which  will  take  over  the  function  of  beam 
control  and  allow  the  focus  anode  to  be  run  at  or  near  cathode  potential. 
Problem  (1)  is  not  being  pursued  at  this  time  because  the  flexibility  of 
the  new  geometry  should  allow  small  adjustments  to  be  made  on  the  focus 
anode  potential,  in  order  to  satisfy  the  Pierce  field  characteristics. 

Note  that  the  intent  of  the  gridded  anode  is  to  shield  the  Pierce  field 
region  from  the  acceleration  region  of  the  gun,  to  prevent  the  lensing 
effects  seen  in  Figure  6-1.  The  beam  control  function  is  a necessary 
side  effect  of  the  primary  shielding  function,  and  is  predicted  to  require 
potentials  of  100-200  volts  relative  to  cathode. 
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Figure  6-1.  Part  of  the  trajectory  solution  for  the  rectangular  beam 
gun  used  in  the  3662  Amplifier. 


Figure  6-2.  The  Pierce  Geometry  Necessary  to  Draw  a Uniform 
Beam  Off  the  Cathode 
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7.0  PROGRAM  FOR  THE  NEXT  PERIOD 

During  the  next  four  month  period,  effort  will  be  directed  at  developing 
a method  of  producing  election  guns  with  improved  cross-sectional  pro- 
files, by  means  of  increasing  the  accuracy  of  the  piece  parts  and  the 
fixturing  used  for  assembly. 

A grid  will  be  added  to  the  first  anode  and  the  gridded  gun  will  be  con- 
structed and  tested.  The  gridded  design  should  provide  a more  well- 
shaped beam  because  emission  can  take  place  over  the  entire  cathode 
surface  at  all  levels  of  cathode  current. 

Another  modification  which  will  be  done  to  the  existing  gun  is  to  increase 
the  length  of  the  anodes  and  focus  electrode,  leaving  the  cathode  at 
0.700  inches.  This  will  improve  the  beam  shape  near  the  ends  of  the 
beam. 

The  slit  diode  array  having  0.002  inch  wide  slits  will  be  fabricated  and 
calibrated,  and  will  then  be  used  to  test  EBS  sheet  beam  guns.  The 
results  of  the  beam  cross-section  analysis  will  be  correlated  with  mea- 
sured intermodulation  performance  of  the  same  electron  gun. 
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GUNOOOIO 
GU.N0002  0 

c 

GUN00030 

c 

GUN00040 

c 

RCCT ANGULAR  GUN  ANALYSIS  PACKAGE 

GUN00050 

c 

GUN0006  0 

c 
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GUN0007  0 

c 

GUN0008  0 

c 

GUN00090 

c 

MAIN  ORIVEP  ROUTINE,  PART  I - EQU IPOTENTI AL  SOLUTIONS 

gunooiog 

c 

GUN001 1 0 

IMPLICIT  REAL *8  (A-H,0-Z> 

GUN00120 

REAL*8  V(41,401) 

GUN001 30 

REAL*3  CUN0I9) 

GUN00L40 

REAL*8  XI 200) , Y ( 200 ) , UX ( 200 ) ,UY( 200) 

GUN00150 

RFAL*8  TITLE ( A ) 

GJN00160 

INTEGER  IC0NDI9) 

GUN001 70 

COMMON  /A/  V 

GUN00180 

CCM.MCN  /R/  XM  IN, XMAX,DXO, 

GUN DO  19  0 

* YMIN,YMAX,DYO, 

GUN00700 

* M, N ,NXY 

GUN0021 0 

COMMON  / C/  CCNC , ICUND, NCOND, CM IN  » C MAX ,COFF 

GUN0022  0 

LCMMON  /□/  X,Y,UX,UY 

GUN00230 

C CM MON  /E/  TITLE 

GUN00240 

DATA  Bw  / .2/ 

GUN00250 

OATA  NEQP  /3 1/ 

GUN00280 

DATA  YEOP  /4.SC-3/ 

GUN00270 

DATA  RELERR  /l. 0-4/ 

GUN0028  0 

OATA  ITER  /20/ 

GUN00290 

DATA  X l , X2 , Y1 , Y2  /O., 4. 0—3,0. ,6.0-3/ 

GUN0030  0 

DATA  ENGMCT  /39.37/ 

GUN003I 0 

c 

GUNJ0320 

c 

RE  AO  IN  INITIAL  INFORMATION 

GUN00330 

c 

GUN00340 

CALL  INIT 

GUN00350 

c 

GUN00360 

c 

READ  IN  VOLTAGFS  FROM  PREVIOUS  RUN 

GUN00370 

c 

GUN00380 

READI8)  V 

GUN0039  0 

c 

GUN00400 

c 

READ  IN  CGNCUCTOR  CONFIGURATION 

GUN004 1 0 

c 

GUN00420 

CALL  MFSH 

GUNJ0430 

c 

GUN0044  0 

c 

RFLAX  VOLTAGE  MATR  I X 

GUNJ0450 

c 

GUN0046  0 

ERR MAX= RELERR* (CMAX-CM  IN) 

GUN904  7 0 

CALL  RELAX  ( 1 , 1 , ERROR , ERRD AX  ) 

GUN00480 

IF  (ERROR. LT.EPRMAX)  GO  TO  A 

GJN00490 

KSTCP=  l 

GUN00900 

1 

K ST  EP=2  *KST  E P 

GUN00810 

IF  I M0D(M-1»KSTEP ) .NE.O)  GO  Tn  2 

GUN00920 

IF  (MOO(N-1,KSTFP).NE.O)  GO  TO  ? 

GUN00530 

GC  TO  1 

GUN00940 

2 

KSTEP=KSTEP/2 

GUNJ0S50 
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CALL  RELAX  ( I TER, KSTEP , ERROR , ERR  MAX) 
IF  (ERRHR.LT.  ERRMAXI  GU  TO  3 
CALL  RELAX  ( I T ER , l , ERROR , CRRM A X » 

RECORD  V MATRIX  FOR  NEXT  USAGE 


REWIND  8 
W R l T E ( 8 ) V 
WR I TE ( l .100) 


FRROR 


PLOT  CONTOURS 
CALL  CCNPLT  (fiW.NEQP) 

DETERMINE  f GUIPOTENTIAL  SURFACE_IN  VICINITY  OF  Y= YEQP 

VO=PINT (Xl,Yl+.75*( Y2-Y1) )-COFF 
CALL  EQPTL  ( VO , XI , X2 , Y l , Y2 ) 

WRITE  OUT  THEN  CONVERT  TO  INCHES 

WRITE(l.lOl)  VO 

OH  10  1=1, NXY 

WR I TE ( 1 . 102 ) X(I),Y(I) 

XU)  = X<  I ) * ENGMET 
Y ( I ) = Y(  I ) AENGMET 
CONTINUE 

GENERATE  INPUT  FILE  FOR  XMGUN 

WRITE ( 2,103)  TITLE! 3), VO, VC 
X ( 1 ) = X (2) 

Y ( l ) = . 2 1 2 *> 

NCARDS=lMNXY-l)/7 

WRITE (2, 10 A)  NCAR0S,C0ND(2) , CONO ( 3) 

DC  20  1=1, NXY, 7 

JF  = MIN0(8,NXY*-l-I ) 


WRITE (2, l 05) 
CONTINUE 
WRITF( 2, 106) 
STOP 

FORMAT  (///• 
FORMAT  <///• 


( Y ( NX  Y*  2-  I- J I , X (NXY  2-  I- J ) , J=l,JF) 


FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 


(///•  MAXIMUM  CHANGE  IN  MESH  ON  LAST  ITERATION  = ',012. A) 
(///'  COORDINATES  IN  METERS  OF  EQUI PO TENT  I AL  LINE  », 

• VO  = • .D12.A//6X, • X', 1LX, • Y ' / ) 

(20 12. A) 

( A8,2X,2F10.0, • .10  622  100  60  1000 100 1 1 00 • ,6X 

•6. •/*. 21251000. .0115. 015  .0115.050  .0115.100  .0115') 

( * 3* , 9X , 1 2, ' 2* ,6X,2( • l'F9.0)  ) 

(8 (2F5. A) ) 

C.O  .200  .0  .017  .01A9.017  .A2911.017* 

/• .AIA21.081.10A1.125A.10A1.0  807.119  . 080 7 . AA391 . 081  • 
/3t*.0*,8XI,*l.  .0  -.009875*) 


GUN00560 
GUN00570 
GUN00580 
GUN0059  0 
GUN0060  0 
GUN006 1 0 
GUN0062  0 
GUN00630 
GUN006A0 
GUN00650 
GUN00660 
GUN0067  0 
GUN00680 
GUN00690 
GUN00700 
GUN907 10 
GUN00720 
GUN00730 
GUN007A0 
GUN0075C 
GUN00760 
GUN00770 
GUN00780 
GUN00790 
GUN00800 
GUN00810 
GUN0082  0 
GUN00830 
GUN008A0 
GUN00850 
GUN0086  0 
GUN0087  0 
GUN00880 
GUN0089  C 
GUN00900 
GUN009 1 0 
GUN0092  0 
GUN00930 
GUN009A0 
GUN0095  0 
GUN0096  0 
GUN00970 
GUN00980 
GUN00990 
GUN01000 
, GUNO 1010 
GUN01020 
GUN01030 
GUN010A0 
GUNO 1050 
GUN01060 
GUN01070 
GUNO 108  0 
INI  000 1 0 
INI  00020 


> 
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SUBROUTINE  INIT 
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CATA  INITIALIZATION  ROUTINE 

IMPLICIT  REAL*£  <A-H,C-Z) 

REAL*8  V(41,401) 

RE AL*fl  CONOm 
R C A L *3  TITLEI4) 

INTEGCR  ICOND(S) 

INTEGER  NCHAR I S ) » L I NE ( 120) . SPACE 

DATA  NCHAR  / 1H  1 , 1H2,  1H 3,  1 H4 , 1 H5 , IH6,  1H7, 1H8 1 IH9/ 

DATA  SPACE  /LH  / 

CCMMON  /A/  V 

CCMMCN  /B/  XMIN, XMAXtOXO, 

* YMIN, YMAX.DYO, 

» M» N»NXY 

COMMON  /C/  CCNC»  ICONO*  NCOND. CM  IN . C.MAX.  COFF 
CCMMON  /E/  TITLE 
OAT  A CMINO  /I.C4/ 

CATA  DECK  SETUP  ON  LOGICAL  RECORD  3 (ALL  UNITS  MKS ) 


CARO 


VARIABLE 


T ITLEd-2) 


T I TLE ( 3 -4 ) 


DXO.DYO 


FORMAT 


DESCRIPTION 


IDrNTIF ICATION 


TITLE 


V MESH  SIZE 


2D  I 2. 2 V MESH  INCREMENTS  ( M ) 
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ICCND.COND  I1,3X,  CONDUCTOR  COOES  - ICCNJ  IS 

r 6.0  THE  CHARACTER  USED  IN  THE  MESH 
TC  REPRESENT  POTENTIAL  CON!) 

(9  OP  LESS  MAY  3E  SPECIFIED) 

(END  OF  RECORD  CARD ) 


R EAD ( 3 . 100)  (TITL£<  I),  I»1,4),M,N 
FEA0(3,101>  DXO.DYO 

ESTABLISH  BOUNDARIES 

X M I N=0 . 

XMAX  = XMIN*OXO*CFLCAT(,M-l ) 

YMIN*0. 

YVAX=YMIN*OYO*CFLCAT(N-I) 


I N 100030 
INI0004C 
IN  100050 
INI00060 
I N 1 0007  C 
INI00080 
INI  00090 
I N 1 00 100 
INI  00 1 1 0 
I N 1 00 12  0 
INI  00130 
IN  I JO  1 40 
INI  00 1 50 
IN  I 00160 
INI  00170 
I N 1 00 18  0 
I N I 00 19  0 
1NI00200 
INI  00210 
IN 100220 
INI  0023  0 
I M 0024  0 
INI  00250 
I N I 0026  0 
I N I 002  7 0 
INI00280 
INI  00290 
INI00300 
INI  003 10 
IN I 00  320 
I N I 0 0 33  0 
INI 00340 
INI00350 
INI  00360 
I N 100  370 
I N 1 00  3ttO 
I N 1 00  390 
INI  0040 0 
INI  004 1 C 
INI  0042  0 
INI  00430 
INI  0044  0 
1 N 10045  C 
INI00460 
I N 1 004  7 0 
INI004eC 
INI  0049  0 
I N 1 00500 
I N 1 005  I 0 
INI  0052  C 
I N 1 0053  0 
INIQ0540 
INI00550 
I N 10056  0 
INI  005  7 0 
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c INI00580 

C WRITE  OUT  GEOMETRICAL  AND  INITIAL  CONDITION  INFORMATION  INI00590 

C _ _ _ _ INI00600 

WRITE! 1,102)  (TITLE! II  * I*i#4)»M»N»  INI00610 

* XMINtXMAXf DXOf  INI00620 

*  YM I N »_YMA  X , D YO  _ INI00630 

C " INI00640 

C READ  CONDUCTOR  CODING  AND  POTENTIALS  INI00650 

C _ _ IN  1 00660 

DO  10  1*1*9  INI00670 

NCOND* I INI00680 

READ ( 3 » 103# END* 1 1 ) ICOND! 1 1 ,COND ( I ) _ __  INI00690 

10  CONTINUE ~ """  - INI00700 

C INI00710 

11  NCOND=NCOND— 1 __  _ _ INI00720 

C - - - ~ INI00730 

C CALCULATE  MAX,  MIN,  AND  OFFSET  INI00740 

C __  ___  INI00750 

CMIN=COND(l)  " INI00760 

CMAX=CONO ! 1 ) INI00770 

C _ INI00780 

DO  20  I*  l , NCOND  " "*  “ * INI00790 

IF  (CMIN.GT.COND! I) ) CMIN*COND!I>  INI00800 

IF  (CMAX.LT.CONC!  !))  CMAX=CONO! I)  __  INI00810 

20"  ‘ CONTINUE"  " ‘ INIG0820 

CCF F*CMINO— CMIN  INI00830 

RETURN  IN  100840 

C "*  " ~ - INI00850 

100  FORMAT  (2A8/2A8/2I5)  INI00860 

101  FORMAT  (3012.41  INI00870 

102  FCRMAT  I/7/m<JX,TDENT  ' " "»,2A87/10X,*TirLF  ' ~',2A8///  INI00880 

* • MESH  SIZE  IS  ',13,'  * * , 13,  ' POINTS*//  INI00890 

* « XM IN  = * ,D10.4,10X,*XMAX  = * , 010. 4, 10X, • XINC  * ',010.4,  INI00900 

* "METERS '/T  — " TNr009l0 

* • YMIN  * «;D10.4,10X,»YMAX  * • , DIO. 4, l OX, ' YINC  * *,D10.4,  IN100920 

* • METERS'//)  INI  00930 

103~"  FORMAT  (T1,3X7F6.0T“  " 1NI00940 

END  INI00950 

C MESOOOIO 

C MES00020 

C MES00030 

SUBROUTINE  MESH  MES00040 

C - - - - . ..  MES00050 

C RECTANGULAR  GUN  ANALYSIS  PACKAGE  MES00060 

C MES00070 

C JOHN  8 RE mC  1/31/79  ‘ - MES00080 

C MES00090 

C THIS  ROUTINE  SETS  UP  A BRAND  NEW  MESH  OF  CONDUCTORS  IN  THE  MES00100 

C POTENTIAL  MATRIX  V.  ‘ ADVANTAGE  IS  TAKEN  OF  ANY  PREVIOUS  MES00110 

C SOLUTION  OF  A SIMILAR  TYPE  BY  CHANGING  ONLY  THCSE  MESH  MES00120 

C VALUES  WHERE  A NEW  CONDUCTOR  IS  LOCATED.  MES00130 

L - - ---  - ...  . MES00140 

IMPLICIT  REAL*8  (A-H.O-Z)  MES00150 

REAL*B  V( 41,401)  MES00160 

«FAL*8  CONDI 9 1 ' " ■'  MES00170 
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noon  non 


INTEGER  I COND ( 9 ) 

INTFGER  NCHAR(S) ,LINE(  120)  , SPACE 

DATA  NCHAR  / IF  1 , 1F2 » IH3 , IH4, IH5, 1 H6, 1H 7 , IH8, IH9/ 

DATA  SPACE  /1H  / 

COMMON  /A/  V 

CCMMCN  /rt/  XM[N,XMAX,DXO, 

* YMIN,YMAX,DYO, 

* M,N,NXY 

CCMMCN  /C/  CQNC » ICOND, NCOND, CM  IN , C MAX ,COFF 
C 

C DATA  DECK  SETUP  ON  LOGICAL  RECORD  4 (CONDUCTOR  GECME  TRY) 

C 

C CARD  VARIABLE  FORMAT  DESCRIPTION 

C 

r.  l-N  M*  A 1 CONDUCTOR  MESH  CCDES 

C 

C LAST  (END  OF  RECORD  CARD) 

c 

c 

C READ  IN  CONDUCTOR  INFORMATION 

C 

00  20  J=  1 » N 

P.FADI4 , 100)  (L INEI I ) , I=l,M) 

DO  10  1=1,  M 

V(I,J)  = OMAXKDABSIVI  I,J)  JtCMINO) 

IF  (L  INEU>.EQ. SPACE)  GO  TO  ID 
DC  11  K=l,9 

KO  = K 

IF  (LINE(I).Eg.NCHAR(KO) ) GO  TO  12 

11  CONTINUE 
GO  TO  30 

12  DO  13  K= l » NCOND 
K 1 = K 

IF(KO.CQ.ICONOCKl) ) GO  TO  14 

13  CONTINUE 
GC  TC  30 

14  VII, J)=  -<C0N0(K1  M-COEF) 
n CCNTINUE 

?)  CONTINUE 
R r T U 0 N 

3D  WRITE!  1,101)  J,  I , ( L INEI  I ) , I = 1 ,M ) 

STOP 

C 

100  FORMAT  ( i 2 0 A 1 ) 

101  FORMAT  (/////•  INPUT  MESH  CRkOR'.lCX,'  LINE  ',15, 

* 10X,'  CHARACTER  • , l *5//l  X , 12  0A1) 

L:NO 
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MES00180 
ME  S00190 
ME  S00200 
MES00210 
ME  S00220 
ME  S00230 
MES00240 
MES0025C 
ME  S0026  0 
ME  S 002  7 0 
MFS0028J 
MES0U290 
•^ES  00300 
ME  SOU  310 
MES00320 
MES00330 
MES00340 
ME  S00350 
MFS00360 
ME S00370 
MCS003R0 
ME  SO 039  0 
ME  S00400 
ME  S004  l 0 
MFSU0420 
ME  S 0043  0 
M6S00440 
MES00450 
MFS00460 
ME  SO 04 7 0 
ME  S0048  0 
MES0049C 
M F S 0 0 9 0 0 
MESU051C 
MES00S20 
MES00S30 
MESOOS43 
MESOOjSO 
MCS0056C 
ME  S005  7 0 
MES00S80 
ME  S9059  0 
MESJOoOO 
MFS00610 
ME  S00o2  0 
MES00630 
MFS0064C 
RCLOOOIO 
REL0002D 
RFL  30030 
REL00040 
PEL  00050 
REL00D60 
REL00070 
RELOO  )80 
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C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

r 

n. 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

r. 

c 

c 


r. 


-NCDF  JISCRTTIZED  LAPLACIAN 

CPEk  ATCK 
M*N  MESH. 

EMPLOYED  TC 
THAT  STARTS  WITH 
THE  FINE  MESH  (I.E.  EVERY 

CF  MESH  GRACES  EMPLOYED  IS  SPECIFIED  BY  K STEP  , 

RESTRICTION  THAT  KSTEP  tiC  AN  INTEGRAL  DIVISOR  CF  BOTH  (M-1I 
AND  ( N — 1 ) . AFTER  EACH  COARSE  RELAXATION  IS  COMPLETE, 

NEXT  FINER  MESH  IS  INTERPOLATED  LINEARLY.  THIS 
PERFORMED  USING  A FOUR  NFAREST  NEIGHBOR 


THIS  ROUT  INF  IMPLEMENTS  THE  4- 

TO  SOLVE  LAPLACE'S  EJUATIGN  WITHIN  A SPECIFIED 
A SUCCESSIVE  CVERRElAXATION  TECHNIULE  IS 
PRCVIDE  FAST  CONVERGENCE,  ALONG  WITH  A SCHEME 
A COARSE  MESH  AND  SUCCESSIVELY  wCRKS  TO 

SINGLE  MESH  PUINTI.  THE  LARGEST 

WITH  THE 


IS 


THE 

I NTLRPCLATICN 
CALCULATION  f MST 
FOR  T 30SF  POINTS  FALLING  AT  A I STANCE  SORT ( 2 ) FROM  THE  COARSE 
MESH,  AND  THEN  FUR  THOSE  POINTS  FALLING  AT  DISTANCE  1.5} 

FROM  THE  COARSE  MESH. 

THE  LAST  FEW  STEPS  FOR  THE  FINEST  GRADE  MESH 
RELAXATION  ARE  GAUSS-SFIOEL  (W=.25l  IN  ORDER 


THr  ROUGH  EDGES  LEFT  OVER  FROM  THE  C VC RRE L A X A T I LN , 


MORt  CCARSF 


WILL 

MESHES 


IT  K S T E P0-  I , OMY  K STE  P=  1 
REPEAT  UNTIL  FRRDp<ERRHAX 
USED  TO  START,  WORKING  TO  THF  FINER  MESH  AND  REPEATING  THIS 
MESH  UNTIL  ERROROTRRMAX 
4l«“t01  MESH  TAKES  ABOUT  4-5  MINUTES 
fc'RRMAX  SET  TO  .Dt  PFKCFJT  CF  THE  MAXIMUM  LESS 
POTENT  IALS. 


FINE 
AS  A 
*1  TH 
MESH 


REL00090 
RFL DOIOO 
RE  L JO  1 1 0 
RFL00I2U 
REL00I30 
RE  L 00 1 40 
RELOO  150 
PEL JOloO 
PEL001 70 
PEL  )0  1 ft  J 
RELOO 191 
REL00200 

p flu 02 1 o 

PEL  0022  0 
R£l002JD 

P EL  0024  0 

PEL  002 5 0 
RtL002oC 
RFL0C270 
PEL  00280 
PEL  0029 0 
RELOO  300 
RELOCUC 

USE  CAUTION  WHFN  CHOUSING  ERR  MAX  RE LOG  l2U 


( KSTFP^l ) 
Tt  SMOOTH 


OUT 


S USED,  AND 
OTHER  W I SE-  , 


THT  PROGRAM 


ART 


TN  AN  MM 


3 70/14  1 
MINIMUM 


St  E CARNAHAN,  LUTHER,  AND  WILKES,  APPLIED  NUMERICAL 
METHODS,  WILEY, 1969,  OR.  FORSYTHE  AND  WASOW,  FIMTE  D IF  TE B- 
ENCT  METHODS  FCP  PDE,  wILEY,  1060. 

ir  THE  MESH  IS  NOT  ENCLOSE')  BY  CONDUCTORS,  PERIODICITY  UF 
THE  SOLUTICN  IS  ASSJMFtJ  AND  T IE  BUUNJARItS  ARE  FOLOFO  CVER 
(I.E.  «IRRCR  IMAGE  SYMMETRY  IS  ASSUMED).  IF  THIS  CANNOT 
BE  TOLERATED,  BE  SURE  TG  SPECIFY  CUNOUCTORS  ALL  THE  *AY 
ARCIN')  THF  MESF. 

IMPLICIT  RE  AL  * 0 (A-H,0-/) 

R F A L * R V(4l,40I) 

COMMON  /A/  V 

COMMON  /:!/  XMINtXMAX,JXO, 


REL00330 
RELOO 340 
RELOO 350 
RFL0D360 
RCL003  70 
cELOOJftO 
RELOO 190 
RELJU400 
RCLD0410 
PEL0J420 
RFL )04  3 0 
PtL0044C 
RELJ0450 
RFL  00460 
PELU0470 
REL004H0 
REL00490 
PLL00400 


* YMIN,YVAX,.DYO, 

* M, N ,NX Y 

MrSH  CCARSC  NFSS  L0QO 

KSTEP=2*KSTFP0 
KSTEP*MAXO( L,KSlEP/2) 

KCF  F = K S TCP 

E-3.141  59/DFLC  AT  ( M I NO  I M/KSTE.n,N/KSTEP)  » 
w=0 .5/ ( l • *DS I N ( T ) ) 

IF  (KSTEPO.EC.I)  W- . 2 5 
1 *'  A X = I TCP/KSI  EP 

ID*  1 


PEL00510 
REL  0052  0 
PEL 00530 
PEL  00540 
RFL0055  0 
REL00560 
REL  0057  0 
REL J0580 
RtL00590 
RtL0060  0 
PFL00610 
PEL00o2  0 
REL JD43D 
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ri  o o ono  nooooonooooo  ono 


r 


tr=v 

JO=l 
JF  = N 

RELAXATION  LOOP 

r)L  30  L=l,LMAX 

L pKCR=0 . 

OC  20  1=10, ir  , KSTFP 

ILPM  ♦KflFF 
I DN= I -KCF  F 

IF  M.EQ.M)  IUP  = IPN 
IF  <I.rQ.ll  I CN= IUP 
DC  10  J= JO, JF , KSTFP 
VO=V< I , J) 

lFIVO.LT.O.)  GC  TC  10 

JCP  = J«-KUFF 
JCN=J-KQFF 

IF  (J.EC.NI  J C P = J 0 N 
IF  ( J . r C • I ) J!~N  = JUP 

VCFF=W*<  DA8S(  V < IDA,  J ) )-V  D*- DAi'SI  V ( lUP  , J ) )-V0 
* ♦•DABS!  V < I , JDN)  )-V0*  )ABS( V(  I ,JUP)  ) — V 0 1 

V ( I , J)  =VO«-VCFF 

FRKCP=  ) 0 A X 1 (CRKCP,DABS( VOFF ) ) 

10  CfNTINjF 
20  CONTINUE 

•»P  I TE  ( 6 , 200 ) L .KSTFP, W, ERROR 

RETURN  IF  ERROR  SATISFIED  WITH  FINEST  MESH  AND  GAUSS-SEIOEL 
HAS  BEEN  PERFORMED 

IF  (KSTFR. Eg. l.ANO. ERROR. LT.FRR.'IAX. AND. W.FQ.  .25)  RETJRN 
IF  ALMGST  COMPLETE , SWITCH  TO  GAUSS-SFIDEL 
IF  IERROR.LT. 10. *ERRMAX. AN3.KSTEP.EQ.il  W=.25 
30  CONTINUE 

CHECK  IF  INTERPOLATION  NECESSARY  FUR  NEXT  FINER  MESH 
CR  IF  COMPLFTEC 

IF  (KSTEPO.EQ.il  RETURN 
Ic  (KSTTP.FQ. 1 I GO  TO  I 

I NT FRPOL AT C IN6CTWFGN  POINTS 

KCFF=KSTEP/2 
I 0=  I ♦•K OFF 
ir=M-Korr 
jo=  i*Korr 
Jf=A-Kl'FF 


REL00640 
REL006SO 
REL00660 
PEL00670 
REL00630 
REL0C690 
REL00700 
REL007IC 
RELOC720 
REL00730 
REL00740 
REL00750 
RELOO  76  0 
PEL00770 
PEL00700 
PFL00790 
REL 00800 
RELOCdlO 
RELOO^D 
REL00630 
RFL00840 
REL00850 
RFL00960 
REL00870 
RFL00880 
REL00890 
RELJ0900 
PEL 009 1 0 
REL0092C 
PCL00930 
PEL00940 
REL0095  0 
REL00960 
REL00970 
RFL00980 
REL00990 
REL01000 
RFL0I010 
PELOID 2D 
RELUI030 
RELOI04U 
REL01050 
P FLO 1060 
REL01070 
PEL01080 
RE  L01090 
RELOIIOJ 
REL01U0 
REL0II20 
PCL01  130 
RELOIIAO 
REL0I150 
RCL01I60 
REL01170 
RELOUBO 
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OC  60  J=JO,JF,KSTEP 
Jl/P  = J«-KCFr 
JDN= J-KflFF 

DC  40  1=10,  IF  ,KSTEP 

IF  ( V ( I, J). LT.O.)  GO  TO  40 

IUP  = H-KCFF 
I CN= l-KCFF 

V<  I ,J)  = .J5*(U4CS{  V(  IDNtJDM  ) +0A.1SI  V(  I UP,  JON)  ) 

* *D4fiS(V<  10N,JUP) )+DABS( V(  IUP,JUP))) 

40  CCNTINUL 
50  CCNTINJE 
l 0= 1+KCFF 

i r= v-koff 
jc=i 


c 

c 

DC  70  0 = J 0 , J F , K S T E P 

JUP=J+KCF F 
Ji:\'  = J-Ki3FF 

IF  (J.Ci.l)  JON= JUP 
IT  (J.FQ.N)  JUP= JON 
OC  60  1=10,  IF.KSTFP 

IF ( VI  I , J) .LT.O.  ) GO  TO  60 
ICP=I tKCFF 
1 0N~ I -K  OFF 

V ( I ,J) =.25*(OAES (V( I ON , J ) ) *QA3S( V( IUP,J) I 

* +oAes<v(  itJONM  oabs<v(  i,jup)  n 

60  CCNTINUF 
70  CONTINUE 
10=  1 
i f=v 

J0=  l+KCFF 
JF=N-KCFF 

OC  GO  J= JO, JF  ,KSTEP 
JU?  = J>-KCrF 
JON=J-KCFF 

OC  30  1=  10, l F ,KSTEP 

If ( V( I , JJ .LT.O. ) GO  TC  30 

IUP=I+KCFF 

ICN=  I-KQFF 

IF  (I.LS.y)  I UP= I ON 
IF  (I.FG.L)  I CN=  IUU 

V( I ,J ) =.25* (JAPS ( V< ION, J I ) «-DAdS( V(  IUP, J) ) 

* ♦DAPS ( V ( I, JON)  ) ♦L>A05<  VC  I,JUP)  ) ) 

50  CCNTINUG 
90  CCNTINUE 
OC  TO  I 

200  FORMAT  <•***• ,21  IC,F  10.4, F 10. ? I 

END 

C 

C 

c 

SUHWCUTTNE  CCNPLT  (nw.NF^P) 


RELOI 190 
REL0120C 
RFL0I2I0 
RFL0I220 
REL01230 
REL01240 
REL0I250 
RELOI260 
REL01270 
RELO  1280 
REL  0 1290 
RFL01300 
PEL0131  C 
REL01320 
RELOI 330 
PFL01340 
REL0135G 
RELO 1360 
RELOI 3 70 
RELOI 3 30 
REL01390 
REL01400 
RFL014I0 
PEL01420 
REL01430 
RELO  144  0 
RELO  1450 
RELO 1 46  J 
PEL  01470 
P ELO 1480 
RELO 1450 
RTLOISOO 
RELO  1 5 1 0 
REL  01 52  0 
REL 01530 
REL0154 J 
REL01550 
RELO 1 56  J 
RELO l 5 7 C 
PEL01530 
RE  L 0 1 590 
RELO  1600 
RELOI  6 10 
RcLO 162  C 
REL 01630 
REL 01640 
RELO 16 5 Q 
REL  0 1 66  0 
RELOI 6 70 
REL 01680 
PEL0169C 
CONOOOIO 
C0N0O02O 
CDN00030 
CONOO  040 
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JOHN  3 RETTIG 
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THIS  ROUTINE  GENERATES  A CONTOUR  EUU I POTENT  I A L PLCT  ON  THE 
L T NE  PRINTER.  CONDUCTORS,  COOED  AS  NEGATIVE  NUMBERS  IN  THE 
MESH , ARE  S Y M 3 C L I ZED  BY  DIGITS  FR  0M  1-9,  ACCORDING  TO  THE 
SCHEME  SET  UP  IN  THE  ARRAYS  ICOND  AND  COND.  OPEN  SPACES 
ARE  REPRESENTED  BY  LETTERS  ACCORDING  TO  THE  VALUES  OF  Bw  AND 
NEOP.  THE  PARAMETER  NEQP  IS  THE  NUMBER  OF  DESIRED  EOUI- 
POTENTIALS  (MAXIMUM  AH),  WHILE  8W  IS  A BANDWIDTH  SPECIFIER. 
THE  WAY  THC  PARAMETER  SW  WORKS  IS  BEST  ILLUSTRATED  BY  AN 
EXAMPLE. 

ASSUME  N E 0 P = 4 , BW=.L0,  AND  10  < V < 40 

A 9 < V < 1 l 

B 19  < V < 21 

C 29  < V < 31 

0 3D  < V < 41 

HOWEVER,  IF  3W=.50  WITH  ALL  ELSE  THE  SAME 

A 7.5  < y < 12.5 

8 17.5  <V<  22.5 

C 27.5  < V < 32.5 

D 37.5  <V  <42.5 

THUS,  3W=0  * I L L ’R  I NT  A BLANK  PAGE,  AND  BW=l  WILL  PRINT  A 

SOLID  PAGE  CF  CHARACTERS. 

THE  PARAMETER  NEOP  MIGHT  Be  ADJUSTED  BY  THE  PRCGRAM  TO 
OBTAIN  AESTHETIC  EQUIPGTENT IAL  SPACINGS. 


IMPLICIT  REALMS  (A-H.O-Z) 

REALMS  V( 41 ,40 1 ) 

RE  A L * B CCNDI9) 

INTEGER  ICOND (5) 

COMMON  /A/  V 

COMMON  /B/  X M I N , XM  A X » 0 XO , 

* Y M I N , Y V A X , D Y D , 

* M , N , N X Y 

COMMON  /C/  COND.  ICOND,  NO-OND,  CM  IN 
DATA  MAX  / 120/ 

I NT  EGER  NCHARI 5) , ACHAR (48) .SPACE 
DATA  NCHAR  / IHl, 1H2, IH3, 1H4, 1H5, 
DATA  AC  MAR  / IHA, IMP , IHC, l HD, 1HF, 
«=  IHJ,  IHK,  IHL,  1HM,  1HN, 

* 1HS, 1HT,IHU. 1HV, LHW, 

* LHA, 1HB , IHC, 1H0,  IHE, 

* IHJ,1HK,IHL,IHM, IHN, 

* 1HS,  1HT,  IHO,  IHV,  1H.J, 
CATA  SPACE  /IH  /,  UNION!  /IH*/ 


, cmax.coff 

, UNI DNT , PR  T l 120) 
I H6,  1 H 7, 1H8,  LH9/ 
IHE , IHG, IHH, 

LHP, 1HU , LHR , 

1 H X , 1 H Y , 1 H Z , 

IHE, IHG, IHH, 

LHP, 1HQ , 1 HR , 

1 H X , LHY, IHZ/ 


CON 00050 
CGN0006  0 
CCN00070 
CCN00080 
CON0009  0 
CONOOIOO 
CONUOUO 
CGN00L20 
CONOOIJO 
CONGO  140 
CONOO 1 5 0 
CCN00I6Q 
CGNOOi 70 
CONOO 180 
CONOO 190 
CONOO  20  0 
CONQ02 1 0 
CGN0022  0 
C0N00230 
CCN00240 
CGN0025  0 
CDN002t>0 
CGN00270 
CUN00280 
CGNO0290 
CONGO  300 
CUN00310 
CONOO  320 
CONOO  330 
CONOO 340 
CuNOQ  15 J 
CCM036C 
CONGO  17  0 
CGN0C380 
CCN0035 0 
CCNG0400 
CHN00410 
CGN0042  C 
CONOO  +3  G 
C0N00440 
CON00450 
CCN00460 
CCN004 re 
C0ND048C 
CGN0049C 
CONOO  50  0 
CCN00510 
CONOO  52  0 
CCN00530 
CON0054  0 
CQN00550 
CON  )05t>  0 
CQN00570 
CON 00580 
CUN00590 
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V I NC= ( C^AX-CVIN  )/DFL0AT ( NEQP- 1 ) 

CQNOOfa JO 

V L I NC=  OLOS 10 ( V INC  ) 

CCNJ061C 

I F X P=  I D l NT  ( 1 • 0 l * V L l NC  ) 

CCN0062  0 

I FRCT=  1 01  NT  ( l.Ci*10.**(VLINC-0rLO4T(  ICXP)  ) ) 

CON0063  0 

IF  tlFPCT  .EQ.7.CR.IFRCT.FQ.9)  1FF.CT  = IFRC  T + l 

C0N0Go40 

V [NC=D FLOAT  ( IFRCT)*10.**IEXP 

CUN00650 

I y IN- 1 DINT  C 1.0C1*CMIN/VINC ) 

CON0066  0 

I *AX=l+i  DINT!  . 599  *C  MAX /VINO 

CCN00670 

NE0P=INAX— IMIN+l 

CCN00660 

wRITE( 1, 100) 

CCN00690 

DC  10  1= i»  NCONC 

CCN00700 

W R l T F ( 1,101)  1C0ND ( I ) * CONOl  I > 

CCN0071 0 

10 

CCNT  IN'JF 

CCN0O72O 

WRITE! 1*102)  HW 

CONGO  730 

DC  20  1 = 1, ME  CP 

CON JO  74 J 

VCP  NT  = DFL  CAT ( I* IN +1-1 ) *VINC 

CCN0C75O 

VL0=VC  FNT-V INC + 3to/2 . 

CONGO  76  0 

VHT=VCFNT +VlNC*3W/2. 

CQNOC 770 

WRITE!  1,103)  ACHAR  ! I > , VL  U , VCENT  , VH I 

CCN00760 

20 

CCNT INOF 

CONGO 7 9C 

WRITE! 1,104) 

CONOG  300 

N0=  PIN ) ( y A X » M ) 

CCNJOBIO 

icrF=o 

CGN00320 

40 

nc  60  j=  i , n 

CbN  30630 

oc  50  10=1, NO 

CCNJ0840 

1 = 10  + 1 OFF 

CCNOO  35 0 

PRT ( 10 J=SPACE 

CQN0086  0 

V0=  V ( 1 , J) 

CGNOOH  7 C 

IT  (VO.LT.O.)  GC  TO  51 

CGN00350 

Vl=  (VO-CCFF  )/'J  INC+BW/2. 

CCNd0Q90 

IVl  = niNT(Vl) 

CONGO')  00 

IF!  CABS  (Vl-')FLCAT  ( W1 1 ) .GT.3W  ) GO  TO  50 

CON  00910 

PRT ( 10 J =ACHAR(  (Vl-I  IIN  + D 

CCN0092  0 

GC  TO  50 

CCN00+30 

51 

OC  52  K= 1 , NCOND 

CGN00940 

K 0=  K 

CON 00950 

IF  (-V0  . EC  . CONC  (Kfi)+CQFF)  GC  TO  53 

CCN009O0 

52 

CCNT  IN'JC 

CCNO 09 7 0 

PRT! 1 0 ) = JN  I ON  T 

C0NG098C 

GC  TO  50 

CCNOO  ) JO 

53 

I V 1= I CCNO (KO ) 

CCNOIOOO 

PRT ( TO ) =NCHAR ( IV1  ) 

CCNO 10 IG 

50 

CCNT INdE 

CGN01O20 

WRITE! 1,105)  (PRT! I ) , 1=1, VO) 

C0N01O30 

60 

CCNTINUE 

CONOl 040 

IF  (M0+ IOFF.GE  .M  ) RETURN 

CCNO l 05 C 

I CF  F= I UFF  +M0 

CUNOlObO 

■M 0=  M I NO  ( M - 1 0 F F * v 0 ) 

CCN01070 

WRITE! 1, 104) 

CCNO 1080 

GC  TO  40 

CGN0109C 

100 

F CR NAT  </////•  *****  CCNOOCTOR  CODES  *****•// 

CGNOl 100 

* l OX,'  c.cnc  VOLTS'/  1 

CONOl l l 0 

101 

FCRNAT  (12X,I1,F12.2) 

CGNOl 12  0 

102 

FERMAT  (/////'  PCT T NT  1 AL  CODES  #*****// 

CONOl 130 

♦ l’X,'  WINDOW  = SFfc.4// 

CGNOl 140 

9H1 
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* ' CODE  L On  MTAN  HIGH'/) 

l) 3 FORMAT! 2X, A1,4X, 3F12. 2 ) 

104  FORMAT  (IH1/1HC) 

105  r OR  MAT  ( 1 X « 120A  1 ) 
h NO 


SUORQUT  I NE  EGPTL  ( VO , X 1 , X 2 , Y 1 , Y2 ) 

C 

C RFC  T ANGULAR  GUN  ANALYSIS  PACKAGE 

C 

C JOHN  0 RETTIG  1/31/79 

C 

C 

C THIS  PPOGR  A w LOCATES  AND  PROVIDES  COORDINATES  F C R THE 

C EQlJ  I POTENTIAL  L I Nr  UF  VALUF  VO,  IN  THE  V MESH.  A LINEAR 

C INTERPOLATION  IS  PERFORMED  GET  I*  OE  N MESH  POINTS.  iT  IS 

C ASSUMED  THAT  TFF  VOLT  A GF  MATRIX  V IS  STRICTLY  MCNOTUNl L 

C IN  THF  Y DIRECTION,  WITHIN  REGION  A3CUT  THE  ECU  I PU  TE  \ T I AL 

C LINT. 

C 

IMPLICIT  REAL*1’  <A-H,0-Z) 

RFAL*8  V ( 4 1 , 40  1 ) 

REALMS  X( 2D0) , Y( 200 ) , UX ( 200 ) , U Y( 2 DO ) 

COMMON  /A/  V 

COMMON  /8/  X M I N , X M A X , ; ) X 0 , 

* YM  I N.YVAX, OYJ, 

* M , N , N X Y 
REALMS  CDN0I9) 

INTEGER  ICONO(S) 

CCMMCN  /C/  CO NO, I COND , NCCiNO  , CM  IN , C MAX , CDF  F 
COMMON  /J/  X , Y , LX , UY 

vi=vo+cc;rr 

I l = M A X 0 ( 1 , l + IOINT  < ( Xl-XMIN)  /OXOJ  ) 

I2=MIN0(Mf  1 + 10  INTI  (X2-X.M  IN) /OXOI  ) 

J 1=  MAX  3(1 , 1 + 10  INK  ( Yl-YMIN  J/OYO)  ) 

J2=MIN!)(N,  l MO  INK  ( Y 2—  YM  I N ) / ,D  Y D ) ) 

NXY  = 0 

DC  10  1=11,12 

JON= J 1 
J U P = J 2 

V0N  = 0A3S( V I I , J ON ) )-V l 
VUP=OA«S( V( I , JLP ) )- V 1 
IF  (VON*VUP)  1 , l C , 10 
1 JM|)=  ( JDN+ JUP  ) / 2 

V MO  = D A B S ( V ( I , JMC ) >-Vl 

IF  (JJN.EQ.JMDI  GC  TO  4 
IF  ( V 0 N * V M D ) 2 , 1C , J 

? JUP  = J MO 

V U P = V M 0 
GC  TO  1 

3 U 0 N = J M D 
VDN=VMO 
GC  TD  1 


CCNOl ISO 
C0N01 16C 
CCN01 1T0 
C0N0118J 
C0N01190 
EQPOOO 1 0 
EQP00020 
EGP0003  0 
EGP0G040 
EQP00050 
EQP0006  0 
EwPOOD70 
EJP00080 
EQP00090 
E G ° ) 0 1 0 0 
EGP001 10 
ECP00120 
EQPD0130 
ECP0014 J 
EGPOO  150 
ECP00160 
FGPOO  170 
E G ° D 0 1 8 0 
EGP00190 
FQ POO  20  0 
EQP00210 
ECPOC22C 
EGPJ0230 
FQP00240 
F0PD0250 
EGP00260 
EGPOO 270 
EQP00280 
EGP00290 
EGPOO  30  .3 
EOP0031C 
E 0 P 0 0 3 2 0 
EQP00330 
80*00340 
EGP003S0 
Ew POO  36  0 
EGP00370 
EGP00380 
EQP00390 
ECP00400 
EGP00410 
ECP00420 
EUP00430 
EQP00440 
EGP00450 
EQP00460 
EGP0047C 
FlP00480 
EGP00490 
EGP00500 
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nnononrsnoono  o o 


■ 


4 N X Y = N X ■»  l 
OJ-O. 

IF  r -1N.FG.VUP>  GO  TO  5 
0J~\ ON/  < VUP-VON  ) 

5 X '.N>  » f = XMIN*DFLCAT  ( 1-1  )*0X0 

Y i NXY  ; = YMIN*OYC* (DELCAT( JM0-1  )-0J  ) 

10  CONTINUE 
If.  TURN 
EN  ) 


FUNCTION  PINT  ( X , Y ) 

RECTANGULAR  GUN  ANALYSIS  PACKAGE 
JOHN  .3  RETTIG  1/31/74 


PINT  INTERPOLATES  THE:  V MATkIX  IN  RE  TWEEN  THE  MESH  POINTS 
USING  A 4 POINT  LINEAR  INTERPOLATION.  IT  ANY  BOUNDARY  IS 
EXCRETED,  MIRROR  IMAGE  SYMMETRY  IS  ASSLMEO  A3CLT  THAT 
BOUNDARY  ANC  TEF  INTERPOLAT  ION  PERT  OHM Ei)  AS  IE  THE  POINT 
FELL  IN  THE  MESH. 


IMPLICIT  RFAL*P  (A-H,G-/> 

R E A L * 8 7(41,401) 

COMMON  /A/  V 

COMMON  /B/  XMIN,XMAX,OXO, 

* YMIN, YMAX.OYO, 

* M » N » N X Y 
X INT=X 

Y INT  = Y 

IF  (X.LT.XMINI  X INT=2.*XMIN-X 

IF  (Y.LT.YMINJ  Y INT— 2.*YM IN-Y 

IF  (X.GT.XMAX)  XINT=2.*XMAX-X 

IF  (Y.GT.YMAX)  Y I NT  = 2 . * YM  AX- Y 

<0= EEL  OAT (M-I I *(  X INT-XMIN) /( XMAX-XMIN) 

YC- HFL  0 AT ( N— l ) * ( Y INT-YMIN) / ( Y v AX- Y MIN) 

I X=  1 + IDINT (XU  I 

IY=  1HDINTIY0) 

FX=XO— I3FL0AT ( IX-I ) 
r Y=  YO— OTLO AT ( I Y— 1 I 
VO=CABS (7( IX,  IY)  ) 

V 1 = 0 AOS (V ( IX+l.IYJ) 

V2=DABS  ( V ( IX,  IY*l>) 
V3=0ABS(V(IX+l,lY*l)J 
l)Vl  = VO+  ( Vi-VO  ) *F  X 
0V2=V2+ ( V 3—  V 2 ) * F X 
PINT=DV1>( DV2—HV 1 ) *F  Y 
RETURN 

END 


ECP00510 
EGP0OS2  0 
FJP00530 
EQP00540 
ECP0085  0 
FQP0056  G 
EQP00S7C 
FQP00580 
EQP00500 
PIN00010 
PIN00020 
PIN00030 
PIN00040 
P I NOOO**  0 
PINOO  16  0 
PIN0007C 
PIN0008G 
PINGOOOO 
PIN  30100 
PINOOllO 
PIN 3 0120 
PINOO  130 
PI N00143 
PIN00150 
PIN0016C 
PINOO I 70 
PINOOIBO 
P INTO  190 
PIN  10200 
PIN00210 
PIN 00220 
PIN 00230 
P IN00240 
PIN00250 
PIN00260 
PIN00270 
PIN00280 
P I N0C29  C 
P I NO  0 30  0 
PINOO 310 
PINOO 320 
PIN00330 
PINOO  140 
PI N00350 
PIN00360 
PIN  303  7 0 
PIN00380 
PIN00390 
P1N00400 
P1N00410 
P I NO 042  0 
PIN00433 
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on  noonnnoo  non  nnnonnnoo 


GUN00010 

GUN00020 

RECTANGULAR  GUN  ANALYSIS  PACKAGE 

G JN00030 
GUN00040 

JOHN  B RETT  I G 1/3  1/7*9 

GUN00050 

GUN0006C 

MAIN  DRIVER  ROUTINE,  PART  III-  TRAJECTORY 

CALCULATIONS 

GUN0007  0 
GUN0008  C 

IMPLICIT  REAL*6  <A-H,C!-Z) 

GUN00390 

GUNOOIOO 

REALMS  V(41,40i) 

GUN001 10 

R E AL *8  C0N0(9) 

GUN00120 

REAL *8  XI 200) , YI200) , UX ( 200 ) , UY( 200 ) 

GUN00130 

INTEGER  ICCNDIS) 

GUN30140 

COMMON  /A/  V 

GUNOOiSO 

COMMON  /8/  XMIN,XMAX,0X0» 

GUN 00160 

YM  I NtYMAX,nY0, 

GUN001 70 

M, N ,NXY 

GUNOOISO 

COMMON  /C/  CCNC, ICONO, NCOND »CMIN»CMAX,COFF 

GUN00190 

COMMON  /0 / X»Y  »UX»UY 

GUN00200 

COMMON  /£/  TI TLEl,TITLE2»T  ITLF3»TITLE4 

GUN002 I 0 

CALL  IN  IT 

GUN0022  0 

REA0I8)  V 

GUN00230 

CALL  TRAJ 

GUN0Q24  0 

STOP 

GUN00250 

END 

GUN0026  0 

- 

IN  1 000  L 0 
I N 1 00020 
IN  1000 30 
INI  0004  0 

SUBROUTINE  INIT 

RECTANGULAR  GUN  ANALYSIS  PACKAGE 

I N 1 0 0 0 5 0 
I N 1 0006  0 

JOHN  3 RETTIG  1/31/79 

IN  100070 
IN  100080 

DATA  INI TI ALIZATION  ROUTINE 

IN  100090 
INI  00 100 
IN  1 0011 C 

IMPLICIT  REALMS  (A-H,0-Z> 

INI00120 
INI001 30 

R 0 A L * b V(4l,401) 

IN  100  140 

REAL*8  CO NO (9 ) 

INI JC150 

R F A L * 8 TITLEI4) 

INI  00 160 

INTEGER  ICONO I 9 ) 

INI 00  170 

INTEGER  NCHAP  ( S ) , L I NE I 1 20 ) , SP ACE 

INI00180 

DATA  NtHAR  /IF  1, lh2, 1H3, 1H4, 1E5, 1H6, 1H 7, IH8, 

1H9/ 

INI00I90 

DATA  SPACE  /IN  / 

INI  00200 

CCMMCN  /A/  V 

I N 1002 l 0 

COMMON  /rt/  XMIN, XMAX,DXO, 

I N 100220 

YMIN,YVAX,DYO, 

INI00230 

M , N , NX Y 

INI0C240 

COMMON  /C/  CllNC, ICONO, NCONO,CM IN, CMAX  ,COFF 

I N 1 0025  J 

COMMON  /E/  TITIC 

I N 1 0026  0 

DATA  CMINO  /I .04/ 

INIC0270 

OAT  A DECK  SETUP  ON  LOGICAL  RECORD  3 (ALL 

UNITS  MKS) 

IN  100280 
INI  00290 

II- 1 


ooo  o onn  o o r>  non  oooonononnocionnonon 


CARD 

VAR  TABLE 

FORMAT 

1 

T I TL  E ( 1-2 1 

2 A 8 

2 

TITLE (3-4) 

2 A 8 

3 

M,  N 

215 

4 

DXC.DYO 

2C12.2 

5-13 

ICCNO.COND 

I l , 3X, 

F6.0 


LAST 


RE AO (3, 100)  ( TITLE! I ), 1 = 1, A) ,M,N 
P E A 0 ( 3 » 101)  DX  C, OYC 

ESTABLISH  BCUNCAR  IE  S 


OESCRIPT  ION 
IOENTIF  ICATICN 
TITLE 

V MESH  SIZF 

V MESH  INCREMENTS  (Ml 

CCNOIJCTOR  CODES  - ICGND  IS 
THE  CHARACTER  USED  IN  THE  MESH 
TO  REPRESENT  POTENTIAL  COND 
(9  CR  LESS  MAY  HE  SPECIFIED) 

(END  OF  RECCRD  CARD) 


XMIN=0. 

XMAX  = XMIN  + DXn*CFLEAT (M-l  ) 

YMI,\=n. 

YMAX=YMIN*OYO*CFLCAT{N-l) 

WRITE  OUT  GEOMETRICAL  AND  INITIAL  CONDITION  INFORMATION 

WRITE!  1*10  2)  (TITLC(I)»I=l,4),M»N» 

* XMIN,X‘1AX,  0X0, 

+ YMIN, YMAX.DYO 

READ  CONDUCTOR  CODING  AND  POTENTIALS 

DC  10  1=1,9 

NCCND=  I 

RFADI  3,  103,ENB=ll  ) I CON'D  I I ),CCND(  I ) 

ID  CONTINUE 

11  NCCND  = NCDND- 1 

CAL  0 JL  ATE  MAX,  MIN,  AND  OFFSET 

CM  I N=CONO ( 1) 

CMAX=CGND(1) 

C 

DC  20  1=1, NCUND 

IF  (CM  IN. GT. CONDI  I)  ) CM IN=COND( I ) 

IF  (CMAX.LT.CCND(  I ) ) C MA  X=C0N3(  I ) 

20  CONTINUE 

C.CFF=C  M INO-CMI  N 
RETURN 

II— 2 


I N 1 00  JO 0 
INI00310 
I N I 0032  0 
I N 1 00330 
IN  I (3034  0 
INI  0035  C 
IN  100 36  0 
IN  100  3 70 
INI30J30 
IN  1 30  390 
INI  00400 
IN  1004 1 0 
INI  0042  C 
INI  004  30 
I N I 0044  0 
INI  00450 
INI0C46C 
INI  0047 0 
I N I 0046  0 
IN l 0049 C 
INI0050C 
INI 00510 
I N I 00  520 
IN l 00 5 30 
INI  00540 
IN  1 00  55  0 
I N I 005o  C 
IN  I 905  70 
i N I 00  5 SO 
IN  100590 
IN  10060  0 
IN I006L0 
INI  00620 
I N 1006  3 C 
INI  0064  0 
INI  006 50 
IN  100660 
INI  30670 
INI 30630 
I N I 0069  G 
I N I 00  70  C 
IN  I 00 7 1C 
INI  (JO  72  C 
INI  30  730 
IN  100  740 
INI  JO  750 
IN  10076  0 
INI  00  7 70 
IN  100  780 
I N I 00  790 
I N 1 0080  0 
iNioonio 
INI  0082 0 
INI00330 
IN  I 00340 


— — ■ - ' 


— — 


I 


c 

100 

101 

10  2 


103 


C 

C 

C 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

t 

c 

c 

c 

c 

c 

c 

r. 

c 

c 

c 

c 


rck«u  (2A8/2A8/2I5) 

FCRMAT  13012.4) 

FORMAT  (///// 10X, • I OENT  • ,2A8// 10X, • T I TLE 

• MFSH  SIZF  IS  ’,13,’  * ’,13,'  POINTS’// 

• XMIN  = • ,010.4, 10X,  • XMAX  = • , DIO.  4,  l OX  , • XI NC 

• MF.TEP5*// 

• YM  l N = • , DIO. 4, 1UA, 

• MFTERS’//) 

FCRMAT  (I  l » 3X  » F 6 . C ) 

FND 


10X,  • YMAX  = ',010.4,1: 


'X,  • YINC  = 


• ,2A8/// 

• ,010.4, 
• ,010.4, 


SUBRCUTINF  T R A J 

RCCT  ANGULAR  GUN  ANALYSIS  PACKAGE 
JOHN  fl  RETT  I G 1/31/70 


IF  DESIRED,  PARTICLES  MAY  Pt'  POUNCED  OFF  CF  OCLNDAkICS  ABOUT 
WHICH  THE  FltLC  IS  ASSUMED  TO  BE  MIRROR  IMAGE  SYMMETRICAL. 


CAT  A DECK  SETUP  ON  LOGICAL  RECORD  5 


CARD 

1 


VAR  I ABLE 

X1,Y1,U1,PI 


FORMAT 


DESCRIPTION 


3D  12. 4, 
F 1 2 . 2 


INITIAL  PCSITICN  AND  VELCCITY 
FOR  PARTICLE  1 - XI  AND  Yl  IN 
METERS,  Ul  IN  MFTERS/S , PI  IN 
DEGREES  FRO*  *X  AXIS 


2* 
LAS  T 


(SA^r  FOR  PARTICLE  2,  L TC ) 
(END  Cf  RECCRO  CARJ) 


imrl  ir  I t r r a l * E ( a-h,l-z  ) 

R E A L * 8 X ( 200  ) , Y I 2 0^ ) » U X 1 200 ) , U Y(  200 
CC*«CN  /H/  XMIN, XMAX, CXO, 

♦ Y " I N , Y * A X , f)  Y 0 , 

• M,N,NXY 

CC-<'»''N  /l/  CNF,  I Lf  NO  , N(  0 N , f,  M I N , C yA  X , COF  F 


THIS  PRCGRAM  FINDS  THE  TRAJECTORY  OF  A CHARGED  PARTICLE  IN 
THE  PRESENCE  CF  AN  ELECTRIC  FIELD  CREATED  8Y  CONDUCTORS  CF 
VARIOUS  POTENTIALS.  THERMAL  EFFECTS,  RELATIVITY,  AND  SPACF 
CHARGE  EFFECTS  AkF  NUT  TAKEN  INTO  ACCOUNT.  FIELD  GRADIENTS 
ARF  TAKEN  F Rl M THE  PCTFNTIAL  MATRIX  V (ASSUMED  TC  HAVE  ALREADY 
BEEN  SOL VtO  BY  SOME  MEANS).  THEN,  A DISCRtTIZEC  FOURTH 
CROFR  DISCRETIZED  RUNGF-KUTTA  INTEGRATION  SCHEME  IS  USED  TO 
PIECE  CUT  TFF  PATh  OF  THE  PARTICLE  THROUGH  THE  MATRIX. 


INI  00 35  3 
I N I 00  86  n 
I N I OCd  7 J 
INI0038D 
INI0089C 
I N 1 0 0 7 0 D 
INI  0071 0 
INI  0072  0 
INI00930 
I N 1 00  74  0 
INI  00950 
TRA00010 
TRA00D20 
TRA0J030 
TRA0004  0 
TRA00J50 
TRA00060 
TRAOG J70 
TPA0008 3 
TRA00090 
TRAOOIOO 
TRAOOl IG 
TRA0012C 
TRAOOL 30 
TRAOG 140 
T R A 0 0 1 5 0 
TRA0C16C 
TRA00L70 
TRA00180 
TR AGO  19  0 
TR AO 020  G 
TR AO 02  1 0 
TRA00220 
TRA00230 
TRA00240 
TRA002S0 
TRA0026 0 
TRA0C2  70 
TRA002C  0 
TRAD 0290 
TRA0030C 
TRA0031 0 
TRA0032C 
TRA0CJ30 
TRA00340 
TPA0035G 
TRAOOibC 
TRAOO  370 
TRA0C38D 
TRAOO  390 
TRA00400 
TRA0C41 0 
TR  A0C«,20 
TRA00h3D 
TRAOC-4  ) 


i! 


II- 3 


oooooono  o ooo  non  non  o non 


CCMMON  /J/  XfY»LX»UV 
ft  F AL  *8  F(4)tZ<4) 

FtJlVALfiNCF.  (Xl,Z(l)),(Yl,7(2))#(UXI#Z(3)J,(l,Yl,Z(4n 
INTEGER  PUNGE 
DATA  TMIN  /0./ 

DATA  DT  /2.0-U/ 

DATA  PI  /3.  l4l5926S35D»0/ 

ELECTRON  CHAPGT/^ASSf  HKS 

3 AT  A LTA  /I  .758794601 1/ 

T MAX-T  M IN+2.* ( YYAX-YMIN)/0S2kT ( 2.*ETA*(CMAX-CMINJ » 

1*0 

n i = i ♦ l 

«FAn(r>,  10J,FN0=3C)  XlfYI.UlfPl 

R3  = (PI/ 130. J*P  1 

UXI  = iJl*CCCS(R  J I 

UY1=U1*DS INtRO ) 

x ( i ) = y i 

Y (1) =Y  I 

UX(  1J-UX1 

U Y ( 1 ) = U Y 1 

N XY  = 1 

T = T Y I f. 

CrMf'UTt  STEP 

1 l K = P.UNGE  I^Z.F  ,T,HT  ) 

K = 1 FLAGS  TF AT  3FRIVATIVES  A'F  NCrOE" 

IF  (K.NF.IJ  3C  TF.  12 

CCVPUTC  Gp  A P I T M S IN  POTENTIAL  USING  CENTRAL  DIFFERENCES 

VXH=?INT(X1*.5+CXC, Yi ) 

VXL=3I\T(Xl-.S*:xc,  Yl ) 

7Yb=PINT(Xl,Yl+.5* )YC) 

V Y L = P I NT  (XI  ,Y1-.5*!)YC) 

GRADVX= ( VXH-VXL ) /FXD 
SKADvY - ( YYM-VYL ) / C V D 

COW*  IT  C L)  F R IV  AT  IVES  AN~>  CONTINUE  R UNGF  INTFGR  A T ICN 

- F ( l ) = )X/CT=UX 

- F ( 2 I =DY /!3T- !IY 

- F( 3) =OUX/UT=CTA*GP APVY 

- H 't)=iTUY/L'T  = LTA*GRAi:VY 

F ( l ) = J X 1 

r ( -»)=uy  l 

F <3 ) = £ T A*GP  ADV  X 
* U)rr  T A*CPADVY 


TP.A0045C 
TR A00460 
TRA30470 
TPA00480 
TRA00493 
TRAQG533 
TRA00510 
TRA0052C 
TRA00530 
TRA03540 
TRA33S50 
TRA3056C 
TRA00570 
TRA00580 
TRA00590 
TRAOObOO 
TRA0061C 
TRA03620 
TRA30630 
TRA3364C 
TRA30b5  3 
TRA00660 
TRA00670 
TRAOObPC 
TRAD0690 
TRA 30700 
TRA-30713 
TRA30720 
TRA00730 
TRA00  740 
TRA30750 
TRA Q0760 
TRA007/C 
TRA JC  78  0 
TRA  30790 
TRA00800 
TRA00810 
TRA30S20 
TRA0C33C 
TKA00340 
TRA00350 
TRAPS  360 
TKA0037 3 
T R A 0 0 H R 3 
TRA30R93 
TRA0090C 
T k A009 I 0 
TRA00920 
TRA  3373 J 
TRA3394P 
TRA0075C 
TP AOC  763 
TPA0077C 
Tk A309e  J 
TRAOO  3 


1 1 - 4 


c 

TRA01000 

GC  TO  1 1 

TRA01010 

c 

TRA01020 

12 

NXY  = NXY  + 1 

TRA01030 

C 

TRA01040 

C 

BOUNCE  OFF  CESIREO  BOUNDARIES 

TRA01050 

C 

TR AO  1060 

IF  (Xl.GT.XMIN)  GO  TO  13 

TRA01070 

XI  =2. *XM IN-X 1 

TRA01030 

UXl  =-UXl 

TRA01090 

13 

CONTINUE 

TRA01100 

C 

IF  (Yl.GT.YMIN)  GC  TC  14 

TRA01L10 

c 

YI=2.*YMIN-Y1 

TRA01 12C 

c 

UY l =-UY  1 

TRA01130 

C 14 

CENT  LNUC 

TRA01 1^0 

C 

IF  (XI  .LT.XMAX)  GQ  TO  15 

TRA01  150 

C 

XI  = 2.*X MAX-Xl 

TR AO l 160 

C 

UXl  =-UXl 

TR AO  1 1 70 

C l 5 

CONTINUE 

TRA01 180 

C 

IF  (Yl.LT.YMAX)  GO  TC  It, 

TRA01 190 

c 

Y l=2.*YMAX-Yl 

TRA01230 

c 

UY l = -UY  l 

TR  AO  12  l 0 

C 1 6 

CCNTINUE 

TRA01220 

c 

TRA01230 

c 

RFCORO  POSITIONS 

TRA01240 

c 

TRA31250 

X(NXY )=X1 

TP  AO  126  0 

Y ( NXY ) = Y 1 

TRA0127C 

UX ( NXY ) =UX 1 

TRA01280 

UY ( NXY ) =UY 1 

TP A01290 

c 

TRA01 300 

c 

CHECK  LIMITS 

TRAO 1310 

c 

TRAP  1320 

IF  (XI .LT.XMIN.0R.XI.GT.XMAX.CR.Y1 .L  T.YMIN.0R.Y1 .GT.YMAX 

TRA01330 

* .CR.T.GT.TMAX.CR.NXY.EQ.200)  GCi  TC  20 

TRAO 1 34  C 

GC  TC  11 

TRA01350 

c 

TRA01360 

c 

PRINT  ANO  PLOT  INFORMATION 

TRAO  1 3 7 0 

c 

TP  AO  1 38  0 

2 0 

WRITE! 1,101)  I 

TRA01390 

no  21  K=  1,  NXY 

TRAO 1400 

T=TMI N+OFLOAT ( K-l )*DT 

TR  AO  14  1 C 

U-OSORT  IUX(K)*UX(K)+UY(K)*UY(K)  ) 

TKA01420 

AN3=( l JO./PI )*CATAN2(UY(K ) ,UX(K) > 

TRAO  14 3 0 

WRITE ( 1,102)  T,X(K ),Y (K) ,UtANG 

TRAO  1440 

z\ 

CCNTINUE 

TR A01450 

CALL  T J PLOT 

TR A01460 

GC  TC  LO 

T R AO  14  70 

30 

CALL  PCFF 

TRAOl<»aC 

RETURN 

TRA014YO 

C 

TRA0150C 

100 

FCRMAT  (3012.4,F12.2) 

TRA01510 

101 

FORMAT  (/////10X,*  *****  PARTICLE  *,I2,'  *****'// 

TRA0152  0 

* 6X, *T ( S) • , EX, 'X( w) • , 8X, • Y(M) • , ex, *t(M/S» • »6X» 'ANGLE ' I 

TRA01530 

102 

FORMAT  ( 1X,4C12.3,F 10.2) 

TRA015<,C 

II—  5 


on  no  on  on  on  o onoonnoonnn  non 


FND 


FUNCTION  kJNGE  (N,Y,F,X,H) 

RFC  T ANGULAR  CUN  ANALYSIS  PACKAGE 
JOHN  B RFTTIG  1/31/79 


FOURTH  ORDER  RUNGF-KUTTA  INTEGRATION 

COPIED  FROM  CARNAHAN,  LUTHFK,  £ WILKES,  'APPLIEC  NUMERICAL 
MET  HODS ' , NEW  YORK,  WILCY,  1969,  PP  374-5. 

IMPLICIT  REAL*8  (A-H, C-Z) 

INTEGER  RUNGF 

REAL *8  PH  1(4), SAVFYI4), YIN), FINJ 
DATA  M/0/ 


M = M+- 1 

GC  TO  (1, 2,3, 4, 5), M 

....PASS  1.... 

1 R LNGC  = I 
RETURN 

....PASS  2.... 

2 UC  22  J=l,N 
SAVfcYI JI=Y(J) 

PHI ( J ) = F t J ) 

22  Y (J)=SAVEY(J )+C.5*H#F( J) 

X=X ♦0 . 5 *H 
RUNGE= 1 
RETURN 

....PASS  3.... 

3 DC  33  J=1,N 

PHI (J)=PHI(J)*2.0*r(J) 

33  Y(J)=SAVEY< J)*C.5*H*F( J) 

RUNGE  = l 
RETURN 

....PASS  4.... 

4 DC  44  J=  I ,N 

PHI (J)  = PH I ( J) ♦ 2 . 0*F ( J ) 

44  Y (J ) = S A VE Y ( J ) ♦ H *F ( J ) 

X =X  *0 • 5*H 
RUNGE  = I 
RETURN 

....PASS  5.... 

5 DC  55  J=  l , N 

55  Y(J)  = SAVEY(J)MPHI<  J)*F(J)  )*H/A.O 


TR AO  155  0 
RUN000 1 0 
RUN0002  0 
RUN00030 
RUN00040 
RUN00050 
RUN00060 
RUN00070 
RUN00080 
RUN00090 
RUN00100 
RUN001L0 
RUN00I20 
RUN00130 
RUN00140 
RUN00I50 
RUNOO 160 
RUN0017C 
RUNOO  180 
RUN 00190 
RIJN00200 
RUN002IC 
RUN00220 
RUN00230 
KUN00240 
RUN0025  0 
RUN0026  C 
RUN00270 
PUN0028  0 
RUN00290 
RUNOO  300 
RUNOO  310 
RUN0032  0 
RUN00330 
RUN0034  0 
RUN0C350 
RUN0036  0 
RUN003  7 0 
RUN0038  0 
RUN0039  0 
RUN00400 
RUN004 10 
RUN00420 
RUN00430 
RUN0044  0 
PUN00450 
RUN0046  0 
RUN004  7 0 
RUN00480 
RUN00490 
RUN0050C 
RUN005I0 
RUN00520 
RUN0053C 
RUN0054  0 


1 1- 6 


— 


L. 


I 


! 

I 

I 

: 


■ ! 

4 


*!| 


jl 

■i 


M«0  RUN00S50 

RUNGE=0  RUN00S6C 

RETURN  RUN00570 

C RUN00S80 

EN3  RUN0059C 

C PINOOOIO 

C P I NO 0020 

C PIN00030 

FUNCTION  PINT  (X,Y)  PINJ0040 

C PIN00050 

C RECTANGULAR  GUN  ANALYSIS  PACKAGE  PIN00060 

C PIN00070 

C JOHN  8 PCTTIG  1/31/79  PIN00080 

C PIN00090 

C PINOQIOO 

C PINT  INTERPOLATES  THE  V MATRIX  IN  BETWEEN  THE  MESH  JOINTS  PIN00110 

C USING  A 4 POINT  LINEAR  INTERPOLATION.  IT  ANY  BOUNDARY  IS  P1N00120 

C EXCEEDED,  M IRRCK  IMAGE  SYMMETRY  IS  ASSUMED  AliCLT  THAT  PIN0013D 

C PCJNCARY  AND  TEE  I NTFR PCL A T ION  PERFORMED  AS  IF  THE  POINT  PIN30140 

C FELL  IN  THE  MESH.  PINO0I5O 

C PINOOlbO 

IMPLICIT  REAL*£  (A-H.C-Z)  PIN0017G 

RE  AL *B  V(41,401)  PIN00180 

COMMON  /A/  V PIN0P19C 

COMMON  /B/  XM I N , XMAX ,0X0,  PIN0U20J 

* YMI |\, YMAX,l'YO,  PIN00210 

* M,N,NXY  PINU022U 

X INT=X  P I NOO  ?3 0 

Y I N T =Y  PIN3U240 

IF  (X.LT.XMIN)  X INT  = 2.*XM  IN-X  PIN00250 

IF  (Y.LT.YMIN)  Y INT=2.*YM  IN-Y  PIN002O0 

IF  (X.GT.XMAX ) X INT=2.*XMAX-X  PIN3J270 

IF  I Y.GT. YMAX ) Y INT=2 . * YM AX- Y PIN00280 

XO=DFLC  AT  I M-  I ) *( X 1 N T-X M I N ) / ( XM AX- X MI N ) PI NJ02Y0 

YC  = DFLnAT(N-l  H>1Y INT-YMIN) /( YMAX-YMJ  N)  P I N3C3DP 

I X=  1+ I D INT I XO ) PIN003I0 

I Y=  l+IDINT ( YO  i P I S 0 0 3 2 3 

F X= XO-DFL C AT ( IX-1)  P1N00330 

FY=YO-DELCAT ( I Y-I  ) PIN0034C 

VC=DABS(V(  IX, I Y)  ) P I NOO 350 

Vl  = 0At3S  ( V ( IX*  I , IY  ) ) P I N JO  36  0 

V2=CA8S(V( IX, IYH ) ) PIN0037C 

V3=0A3S(V( IX+l,IY*I) ) PINJC3B0 

OYl-=VO+(Vl-VO)*EX  PINOQj'iO 

DV2=V2+ (V3-V2 ) *FX  P1N004U0 

PINT=DVI*(DV2-DVI )*FY  PINCOhIC 

RETURN  PIN00420 

END  . PIN0043U 

C TJP00010 

C TJP00J2C 

C TJP00D30 

SUBROUTINE  TJPLOT  TJP00040 

C TJP00Ot>0 

C RECTANGULAR  GUN  ANALYSIS  PACKAGE  TJPU0060 

C TJP00070 


1 1— 7 


no  non  non  non  non  ononn 


JOHN  ri  RETTIG  1/31/79 


TRAJECTORY  F*_CT  SETUP  F .’R  THE  ZETA  PL3TTFR 
REAL*1)  V ( 4 1 » 40  1 ) 

RFAL*t  XI 20J) , Y ( 200  ) ,UX(?  JC)  ,UY(  20CI 
REAL*R  XMIN,XvAX,i:XJ,YMri,Y^AX,,)YO,TMIN,rMAX,!')r 
CCMPLEXM6  IDFNT, TITLE 
CCMMCN  /A/  V 

COMMON  /V  XMIN,XMAX,  UXJ, 

* YMIN,YMAX,DY0, 

* M , N , N X Y 

ccmmcn  nt  x , y • lx i uy 

CCMMCN  /E/  IGENTjTITLF 
OAT  A WIDTH  /4./ 

DATA  IFLAG  /0/ 

FIRST  TIME  THROUGH? 

IT  l IF L A J.NE.O)  GC  TO  bO 
I FLA  3=1 

LALl  PLCTF  (10,2) 

CALL  f ACTCR  (2.) 

OEFINF  SCALES 


XL  = w IOTH 

YL=M  DT  H*SNGL ( ( YMAX-YM  IN  ) / ( XMAX-XM  IN  ) ) 
DX=SNCL(XMAX-XYIN)/XL 
0Y  = S.yGL  (YMAX-YV  IN)/YL 
X0=SNGL ( XMIN) 

Y0=SNGL (YMIN) 

TITLE 

CALL  SYMSCL  (YL/3.,XL*l.,.  15,  I"*FNT,0.,  16) 

CALL  SYYOGL  ( YL / 3 . , XL* . 50 , . 1 5 , T I TL E , C . , 16 ) 

SCALE  ALL  FUTURT  PLOTTING  TO  GUN  DIMENSIONS 

CALL  OFFSET  (YC,UY,XO,CX) 

OR  Ai,  CONDUCTORS 
C 

DC  2U  1=1, M 
X 3=  SNGL ( TXO)*FLCAT( I- l ) 

JTEST=0 
DC  10  J=l,N 

IF  (VI  I ,J).LT.C..AND.JTCST.EQ.  0)  GO  TC  11 
IF  ( V ( I ,J).GT.C..ANO.JTCST.EJ.l)  GO  TO  12 
IF  ( J.FQ.N. AND.JTEST.rC. I ) GC  TO  12 
GC  TO  10 
11  JTEST=l 
JSTART=J 


Tjpoooao 
T JP00090 
TJP00100 
TJP00110 
TJP00120 
TJP001 JO 
TJP00140 
TJP00150 
T JP00160 
TJP30170 
TJPOOIRO 
TJPOO  190 
TJP0020C 
TJP0J210 
T JPO  022  0 
TJP00233 
TJP0024C 
TJP0G250 
T JOQ0260 
TJP0027C 
TJP0028C 
TJP3029C 
T JPJ03J0 
TJPOO 310 
TJPOO  J 2 0 
T JPJOJ Jo 
T JP  30  34  0 
TJP00J50 
TJPOO 36 3 
TJPOO  3 7 C 
T JP0G380 
TJP00J90 
TJP0Q400 
T JP0091C 
T J PQ04  2 0 
T JP0043 0 
TJP  )0<*4C 
TJPJC490 
T JP0046U 
TJP 00470 
TJP00480 
T JP  J0'+9  J 
TJPOOSOO 
TJP JOS  10 
TJP 3 0520 
T JPOOS3J 
TJPOO  *>40 
TJP00S50 
TJP00560 
TJP 00570 
TJPOOSt)  3 
TJP00S90 
TJPJ06U0 
TJP0061C 
T J POOF  2 0 
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GC  TC  10 
J T E S T = 0 
JSTCP=J-L 

IF  (J.EC.M  J STGP  = J 

IF( JST ART.EQ. JSTCP)  GG  TO  10 

YO=SNGL(DYO)*ELCAT  ( JST  APT-  l ) 

CALL  PLOTX  (YQ,X3,13I 
YO=SNGL  < )YO)*FLCAT( JSTCP-l  ) 

CALL  PLCTX  C YO  » X 0 , 12) 

CONTINUE 

CCMIN'JE 

OC  40  J=  l » N 
YO=SNGL(DYO)*FLGAT ( J-l  ) 

I TE  ST  = 0 
no  30  I - l » M 

IF  (V( l ,J  I.LT.C..ANO.ITEST.E3.0)  GO 
IF  ( V ( l »J).GT.C..AN3.  ITEST.EO.l)  GO 
IF  ( I .FC.N.ANO. ITEST.EC.l ) GC  TO  32 
GC  TO  30 
I TE  ST=  l 
I ST ART=  I 
GC  TO  30 
I TF  ST  = 0 
I STCP= I - 1 

IF  (I.L«2.M)  l S TCP=  I 

IF(  ISTAPT.EC.  ISTOPI  GC  TG  30 

XO-SNGL  IOXO)*FLCAT  ( I ST  ART- 1) 

CALL  PLCTX  (YO,XO»13) 

XO*SNGL 1 0X0) YFLLAT I I STOP- 1 ) 

CALL  PLCTX  (YOtXO»  12) 

COST  I N(JC 
CONTINUE 

CALL  LINE  (Y,X,NXY, 2,0,12) 

RETURN 

ENO 


T JP00630 
TJP00640 
TJP0Q650 
TJP0066Q 
TJP00670 
TJP00680 
TJPG0690 
T JP00700 
TJP00710 
TJP00720 
T JPC073  0 
TJP00740 
T JP00750 
TJPOO  760 
TJPOO  770 
TJP00780 

TC  31  T JP00790 
TC  32  TJP00800 

TJP008I0 
T JP30820 
TJPOO 830 
TJD00840 
T JP00850 
TJP00860 
T J POOR  7 0 
TJP00380 
TJP00890 
T JP00900 
TJP00910 
T JP00920 
T JP0093C 
T JPJ0940 
T JP00950 
TJP00960 
T JP00970 
T JP00980 
T JP00990 
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Appendix  IV 
KEY  PERSONNEL 


DOUGLAS  B.  CLARK,  Member  of  the  Technical  Staff,  Tube  Division,  Devices 

Group.  Born  December  31,  1946,  Oakland,  California.  Nine  years 
experience.  B.S.,  Electrical  Engineering,  University  of  California 
at  Berkeley,  1969.  Graduate  studies  in  Electrical  Engineering  at 
Loyola  University  of  Los  Angeles. 

Mr.  Clark  is  currently  a Project  Engineer  with  the  Electron 
Bombarded  Semiconductors  Group.  He  is  responsible  for  design  and  develop- 
ment of  state-of-the-art  EBS  devices. 

From  February  1974  to  May  1977,  Mr.  Clark  was  with  the  Solid 
State  West  Division  of  Varian  Associates  in  Palo  Alto,  California.  He  was 
Project  Engineer  on  a variety  of  development  and  production  projects  involving 
Gun  diode  amplifiers,  YTG-tuned  Gunn  Oscillators,  and  bipolar  transistor  amplifiers. 

From  June  1969  to  February  1974,  Mr.  Clark  was  with  the  Hughes 
Aircraft  Company,  Electron  Dynamics  Division,  where  he  was  a Member  of  the 
Technical  Staff.  His  responsibilities  included  that  of  Project  Engineer  on  state- 
of-the-art,  dual  mode,  multi-octave  traveling-wave  tubes.  He  also  was  respon- 
sible for  production  of  high  pulse  power  space  tubes. 


JOHN  B.  RETTIG,  Member  of  the  Technical  Staff,  Tube  Division.  Born 

January  8,  1954,  Toledo,  Ohio.  B.S.E.E.,  1977,  M.S.E.E., 

1978,  Purdue  University. 

Mr.  Rettig  is  presently  working  on  the  development  of  an  EBS  space 
amplifier  for  NASA  Goddard  Research  Center  that  must  meet  stringent  linearity 
and  efficiency  specifications. 

Formerly  he  was  a research  assistant  at  Purdue  University,  engaged 
in  High  Gradient  Magnetic  Separation  studies.  This  involved  experimental  quan- 
tization of  3-dimensional  buildups  of  small  paramagnetic  particles  on  saturated 
ferromagnetic  wires. 

Mr.  Rettig  is  a member  of  Tau  Beta  Pi  and  IEEE. 
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